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Fluctuations in the formation time of ultracold dimers from fermionic atoms 
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We investigate the temporal fluctuations characteristic of the formation of molecular dimers from 
ultracold fermionic atoms via Raman photoassociation. The quantum fluctuations inherent to the 
initial atomic state result in large fluctuations in the passage time from atoms to molecules. As- 
suming degeneracy of kinetic energies of atoms in the strong coupling limit we find that a heuristic 
classical stochastic model yields qualitative agreement with the full quantum treatment in the ini- 
tial stages of the dynamics. We also show that in contrast to the association of atoms into dimers, 
the reverse process of dissociation from a condensate of bosonic dimers exhibits little passage time 
fluctuations. Finally we explore effects due to the non-degeneracy of atomic kinetic energies. 

PACS numbers: 03.75.Lm, 03.75.Ss, 42.50.Ar 
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I. INTRODUCTION 

The coherent formation of ultracold diatomic 
molecules from quantum-degenerate bosonic or fermionic 
atomic gases, via either Feshbach resonances or two- 
photon Raman photoassociation 0], has witnessed 
spectacular developments in recent years and has lead 
to the first realization of molecular condensates Q. Be- 
cause in these experiments the molecular field is initially 
in a vacuum state, it is to be expected that quantum 
fluctuations play a dominant role in the early stages 
of molecule formation. These fluctuations manifest 
themselves in the quantum statistics of the resulting 
molecular field, and also in the time that it takes for the 
number of generated molecules to reach a specific value, 
the so-called passage time statistics. This problem is 
closely related to other situations where quantum (or 
thermal) fluctuations trigger a system to undergo a tran- 
sition away from a dynamically unstable state. One such 
example familiar from quantum optics is superradiance 
(or strictly speaking superfluorescence) , a situation 
where an ensemble of two-level systems initially in their 
excited electronic state and coupled to the electromag- 
netic field vacuum undergoes a transition characterized 
by the emission of an intense light pulse. One important 
difference is that in superfluorescence experiments the 
atoms are normally coupled to a continuum of modes of 
the radiation field, practically leading to an irreversible 
decay to their ground state, while in the problem at 
hand the molecular field is to a good approximation 
single-mode, leading to reversible dynamics. 

One important way to characterize the dynamics of 
molecule formation is by way of the so-called passage 
time, which is defined as the time it takes to produce a 
predetermined number of molecules. The goal of this pa- 
per is to study the passage time statistics resulting from 
the initial quantum fluctuations of the atomic matter- 
wave field. We also compare this situation with the dy- 
namics of dissociation of a molecular condensate into 
fermionic atomic pairs, showing significant qualitative 
differences between the two cases. 



The paper is organized as follows: Section II estab- 
lishes our notation, presents our model, and shows that 
the application of an Anderson mapping leads to the 
description of photoassociation of fermions in terms of 
the inhomogeneously broadened Tavis-Cummings model 
of quantum optics. Section III concentrates on a "ho- 
mogeneously broadened" version of this model that ne- 
glects the spread in fermion energies, an approximation 
shown to be valid for sufficiently small numbers of atoms. 
There we also discuss an approximate stochastic classi- 
cal description 0] that yields a satisfactory qualitative 
agreement with the full quantum results for short enough 
times. This section concludes by comparing the passage 
time statistics associated with photoassociation and the 
reverse process of photodissociation. The results of a nu- 
merical analysis of the full, inhomogeneously broadened 
model are presented in section IV. Finally section V is a 
summary and outlook. 



II. THE MODEL 

In typical experiments that produce molecules via Fes- 
hbach resonance the magnetic field is swept across the 
resonance. In the strong coupling regime kpa > 1 in the 
vicinity of the resonance, where kp is the Fermi wave 
number and a is the s-wave scattering length, the in- 
terpretation of the molecular state is subject to concep- 
tual difficulties stemming from the dressing of the " bare" 
molecular state by atom pairs in the open channel j^. 
Furthermore, since the binding energy of the molecules is 
very small, of the order 10~^^eV, they are larger than the 
interatomic separation. Another consequence of the al- 
most vanishing binding energy is that nearly every time- 
dependent process becomes nonadiabatic with respect to 
the time scale set by the inverse of the binding energy. To 
avoid these difficulties we restrict our considerations in 
this paper to the description of a quantum-degenerate gas 
of fermionic atoms of mass m/ and spin a =t, i, coupled 
coherently to bosonic molecules of mass mi, — 2m f and 
zero momentum via photoassociation rather than Fesh- 
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bach resonance. 

Neglecting collisions between fermions and assuming 
that for short enough times the molecules can be de- 
scribed by a single-mode bosonic field, this system can 
be described by the boson-fermion model Hamiltonian 



|0). We concentrate in the following on times short 
enough that the atomic sample remains essentially un- 
depleted and it is sufficient to consider only fermionic 
levels up to the Fermi surface in Eq. (jSJ. 



k 

+ Tt-xYI {b^Ckfc-ki + ' (1) 

k 

where w , h are molecular bosonic creation and annihila- 
tion operators and cj^^ , Cka are fermionic creation and 
annihilation operators describing atoms of momentum 
hk and spin a. The first and second terms in Eq. 
describe the energy, uju = h^k^/mf, of the atoms, and 
the detuning energy of the molecules respectively, and 
the third term describes the photoassociation of pairs of 
atoms of opposite momentum into molecules. 
Introducing the pseudo-spin operators H 

= ^(clr^fcT + ^-ki^~ki - 1), 

= i^k)^=^lkAv (2) 

which are easily seen to obey the SU(2) commutation 
relations 

[a+,a^,] = 25kk-ijl (3) 
[alai] - ±4fc'<T±, (4) 

where 5kk' is the Kroneckcr delta function, the Hamilto- 
nian becomes, within an unimportant constant [9lll0l|. 

iJ = ^ fiLOk^l + hujbb^b + hxY^ (b^^- + 6(7+) . (5) 

k k 

This Hamiltonian is known in quantum optics as the 
inhomogeneously broadened (or non-degenerate) Tavis- 
Cummings model It describes the coupling of an 

ensemble of two-level atoms to a single-mode electromag- 
netic field. Hence the mapping ^ establishes the formal 
analogy between the problem at hand and Dicke super- 
radiance, with the caveat already mentioned that we are 
dealing with a sing le bosonic mode [S, ilOfc jJ, 111 III 111 ■ 
Instead of real two-level atoms, pairs of fermionic atoms 
are now described as effective two-level systems whose 
ground state corresponds to the absence of a pair, \gk) = 
|Ofc|,0_fc|) and the excited state to a pair of atoms of 
opposite momenta, |efc) = 1-fcx). 

The initial condition of the superradiance problem is 
a sample of inverted two-level atoms. It corresponds in 
the present case to the initial atomic state 

k 

where the product is taken up to the Fermi surface for 
T = 0, while the molecular field is in the vacuum state 



III. DEGENERATE MODEL 

We consider first the simplified situation of a degener- 
ate model in which the inhomogeneous broadening due to 
the spread in atomic kinetic energies is ignored. This is 
justified provided that these energies are small compared 
to the atom- molecule coupling energy, /3 ~ ep/{fix) <C 1, 
where ep is the Fermi energy. This approximation is the 
analog of the homogeneous broadening limit of quantum 
optics, and of the Raman-Nath approximation in atomic 
diffraction. As we will find it is valid only for relatively 
small (~ 10^ — 10'^) particle numbers, but the model ex- 
hibits the essential physics. 

The atom-molecule coupling has been estimated 0, 
m for the case of^^Rbto be ~ 7.6 x 10"^ m^/^s"!, 

so that P = ep/hx « 10^^ N"^/^^. In the last term we 
have related the Fermi energy and volume to the os- 
cillator frequency Who of a spherically symmetric har- 
monic trap via ep ^ N^/'^huJho and V w N^^^a^^^ where 
Oosc = i/Ti/mw/io is the oscillator length. Typical exper- 
iments use traps with Who ~ 100 Hz. In the case of *^Rb 
the trap should contain at most ^ 10^ — 10'^ atom pairs to 
be in this regime, i.e. to have /3 ^ 1. For larger samples, 
it is necessary to account for the inhomogeneous broad- 
ening of the sample, a situation that we consider in the 
next section. For any atom numbers, the characteristic 
time scale is given Tp = 1/xVN ~ 5.8 x 10~^N~^^^s for 
parameters used in this paper. 



A. Quantum description 

Limiting for now our considerations to small atomic 
samples, we approximate all w^'s hy up and introduce 
the collective pseudo-spin operators 

k 

S^-T.-k. (7) 

k 

which again obey SU(2) commutation relations, yielding 
the standard Tavis-Cummings Hamiltonian |llL Il5l| 

H = TiLopS, + hubb'^b + nx{bS+ + 6^5"). (8) 

This Hamiltonian conserves the total spin operator S'^, 
which, by using the pseudo-spin commutation relations, 
can be expressed as 

S+S~ + S,{S,-1), (9) 
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FIG. 1: Short-time dynamics of (hb). From left to right, the 
curves give the linearized solution 12H and the full quantum 
results for iV = 500, iV = 250, and iV = 100, respectively. 

with 

5^1^^) = SiS + 1)\F) = ^(^^ + ?j \F), (10) 
so that S = N/2. Here 

N = b^b + Xl(4T^fcT + cLfcxC_fcx)/2 ^hb + fip (11) 

k 

is the total number of molecules and atomic pairs, which 
is conserved by the Hamiltonian From the definition 
of 5*2 we also have that 

1 iV 1 

S, = - {2hp ^N) = --n,^-in„-hf,), (12) 

hence Sz measures the difference in the numbers of atom 
pairs and molecules. 

Introducing for convenience the joint coherence opera- 
tors 

X = {bS+ + b^S-)/2, 

Jy = {bS+ -PS-)/2i, (13) 

yields the Heisenberg equations of motion 

iib = ~2xJy, (14) 

L - 5Jy (15) 

Jy = -5Jx - X {^S.hb + S+S-^ , (16) 

where 6 ^ LUf, ~ ujp, so that 2xJx + Shh is a constant of 
motion. 

In the following, we confine our discussion to the case 
of 5 = for simplicity. We thus neglect the contribution 
of Jx in Eq. H16() . In order to obtain an analytical so- 
lution valid for short times for (fii,), where () indicates 
the expectation value, and assuming the initial state \F), 
we keep only terms of order fif, on the right in Eq. (|16|l . 
Using Eqs. and (|12|) . we reexpress S^S^ as 

S+S- ^-hl + i2S-l)nb + C^C, (17) 



FIG. 2: Short time dynamics for Anl/{nb). From left to 
right, the four curves give the linearized solution 12211 and the 
full quantum resuhs for N = 500, N = 250, and iV = 100 
respectively. 

where we have introduced for convenience the operator 
(+C ^ - S{S -1). (18) 

Substituting then Eq. H17|) into Eq. H16f) and dropping 
the term proportional to , we have for the early stages 
of molecule formation 

jy^-2xNnh-xC^C- (19) 

Differentiating Eq. (|14|) . taking its expectation value, 
and substituting Eq. (|19|l into the resulting form yields 

fib « {"^Nhb + 2C+r ) (20) 

For the initial state \F) this has the solution 

(Mt)) ~ {{C^C)/n) sinh^ ix^t). (21) 

Similarly, the variance of the molecule number distribu- 
tion is found to be 

Anlit) = (nlit)) - {n,itr 

w ((C+r)/8A^) (cosh (4xViVt) - l) .(22) 

Figure ^ compares the average molecule number (fif,) 
and Fig. H the normalized variance An^/(nb) from Eqs. 
(I21|l and (|22|l respectively, with the full quantum solu- 
tion obtained by direct diagonalization of the Hamilto- 
nian (jHJ. Both approaches agree within 5% until about 
20% of the population of atom pairs has been converted 
into molecules. Note that 

lim Anl/{hb) = 1. (23) 
xvTvt^o 

This is indicative of the fact that for short times the 
molecule field is thermal in character, see for example 
[l0| and references therein. This is further confirmed 
by a comparison of the molecular number statistics to a 
thermal distribution, as illustrated in Fig. |31 



FIG. 3; Molecule number statistics at xV'^i ~ 1-0: solid 
line - full quantum result, triangles - thermal distribution, 
and xVT^t — 2.0: dashed line - full quantum result, circles - 
thermal distribution, for A'^ = 500. 



B. Classical stochastic description 



FIG. 4: Passage time distribution for converting 5% of the 
initial population consisting of only atoms (molecules) into 
molecules (atoms) for A'' — 500. For initially all atoms: solid 
line - full quantum description; dashed line - classical stochas- 
tic model. For initially all molecules: dot-dashed line - full 
quantum result 



Ea. H19|l shows that for the initial state \F) the source 
of the molecular field is the operator C^C • This is the 
non- vanishing part of the second-order moment {S'^S~), 
which is a measure of fluctuations of the atomic field. 
Since the molecular field is initially in a vacuum, its 
growth is therefore triggered solely by these quantum 
fluctuations. 

It is oftentimes possible to simulate the effects of quan- 
tum fluctuations by averaging over a large number of clas- 
sical trajectories triggered by random noise. To imple- 
ment such an approach, we flrst observe that for ^ 1 
the higher-order moments of the pseudo-spin operators 
factorize approximately into a sum of products of 
second-order moments llallSj ^.g., 

{S+S+S-S-) = 2{S+S-){S+S-) + 0{S-^). (24) 

We then proceed by replacing the quantum operator 
by a stochastic classical variable s+ and assuming that its 
fluctuations obey random Gaussian statistics, the prob- 
ability distribution of being given by 



+ ^ - 



1 



/2^A5- 



■ exp 



V2(AS'+)^ 



(25) 



where the variance (AS""*")^ is adjusted to its quantum 
value 



(26) 



for the initial state \F). 

The resulting passage time statistics p{t) can be deter- 
mined numerically by first obtaining classical trajectories 
for (nb(t)) using Eq. (|21|l and assuming that the initial 
values of follow the distribution (|25|l . Taking into 
account that each choice of |s+p within a differential el- 
ement ds~^ maps (fif,(t)) in such a way that it reaches 
a fixed reference value after a uniquely determined 



passage time r within the differential element dr, it fol- 
lows that the probability of a particular value of s"*" is 
equal to the probability of measuring (hb) = nl°^ at that 
time T 

\pis+)\ds+{T,nf) = p{T)dT = \p{s+{T))f{T,nf)\dT, 

(27) 

the differential elements ds+ and dr being related by 



/(r,nr) = |ds+(r,nr)/dr| 



(28) 



By fixing {nb{T)) = nl°^ one can invert Eq. (|21|l which 
upon differentiation with respect to r gives 



(29) 



sinh ^(xV Nt) 



We thus obtain an analytical expression for the passage 
time distribution p{t) as 



Pir) 



' 2rffN 



cosh (xVNt) 
sinh ^ (x 

rof 



Nt) 



exp 



2sinh2(xViVr) 



(30) 



The dashed line in Fig. ^ obtained from Eq. (jJOll, 
shows the distribution of passage times required to pro- 
duce a normalized molecule number rf^^ /N = 0.05 from 
a sample initially containing N = 500 pairs of atomic 
fermions. It should be compared to the solid line, which 
is the result of the full quantum dynamics. The classical 
result reproduces qualitatively the broad and asymmetric 
distribution of the quantum solution. However, it does 
not reproduce well the leading and trailing edges of the 
distribution, which depend on the higher order moments 
of the classical field s+ and are poorly treated by the as- 
sumption of Gaussian noise. In addition the continuous 
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IV. NON-DEGENERATE MODEL 



<n^> 
b 



FIG. 5: EfTective potential for a system with ^ 1. The cir- 
cle (square) corresponds to an initial state with all fermionic 
atoms (molecules). The part of the potential for rii, < is 
unphysical. 



distribution of the classical model fails to properly de- 
scribe the dynamics of the system at very low molecule 
numbers when the discrete nature of molecule number is 
more important. 



C. Photodissociation 

The passage time distribution for the photoassociation 
of fermionic atoms into molecules differs sharply from 
its counterpart for the reverse process of photodissoci- 
ation from a molecular condensate into fermionic atom 
pairs, which is plotted as the dot-dashed line in Fig. ^ 
In contrast to photoassociation, this latter process suf- 
fers significantly reduced fluctuations. One can gain an 
intuitive understanding of this difference by treating the 
short-time molecular population classically, {fit,) —f n^. 
Within this approximation, the Heisenberg equations of 
motion (|14ll- 115|l can be recast in the form of a Newton 
equation jl5(J 



rib 



drih 



(31) 



where the effective potential U{ti{,) is given by 
U{nb) = ^A^^(iV + 3) - 2Nnb - {2N - l)nl + 2ng (32) 

and is cubic in rib, see Fig. O 

In case the system is initially composed solely of 
fermionic atoms, nb{0) ~ 0, the initial state is dynam- 
ically unstable, with fluctuations having a large impact 
in the build-up of Ub- In contrast, when it consists ini- 
tially solely of molecules, Ub = N, the initial state is far 
from the point of unstable equilibrium, and nb simply 
"rolls down" the potential in a manner largely insensi- 
tive to quantum fluctuations. This is a consequence of 
the fact that the bosonic initial state provides a mean 
field that is more amenable to a classical description. 



Although the degenerate model offers the benefit of 
allowing analytical solutions and an intuitive physical in- 
terpretation in the short time limit, we have seen that it 
is only valid for relatively modest values of N. Current 
experiments, however, generally trap ~ 10^ — 10® atoms, 
in which case it is important to properly account for the 
atomic kinetic energies within the Fermi sea. 

From the Hamiltonian jsj, we readily obtain the 
Heisenberg equations of motion: 



drib 

dt 
dt 



k 



-5kil - '^Xnb^k 



-Skit - "^Xfibcrl - Xcrt^k ' 



where we have defined 



Jk 



fk 



2i 



(33) 
(34) 

:) 

(35) 

(36) 
(37) 



and (5fc = — uj^. In Eq. (|35|l we have also made the 
approximation 



Skk'^t<y 



k'"k ' 



(38) 



which is valid for short times when starting from the 
initial state \F). In order to develop a classical model in 
analogy with the degenerate case we define 



I {clck + clfcC^fc) . 



(39) 



Since J2k 



(fib), where the sum runs over momenta 



inside the Fermi sea, the expectation value of this opera- 
tor can be interpreted as that fraction of the initial pair 
of atoms of momenta {^k,k) that has been converted 
into a molecule. Note that 



^k^k +^'ki^k - 1 



with &l\F) = i(i + l)|F) and 



-z 1 'h 



so 



'k "fc 



where we have defined 



^^-a^(<T|-l) 



C^Ck 



1,1 



(o-l)- 



2'2 



(40) 
(41) 

(42) 
(43) 
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FIG. 6: Average molecule number as a function of time for 
various values of the coupling strength. From left to right: 
degenerate analytical, full quantum treatment for /3 = 0.1; 
short-time treatment for /3 = 0.1, short-time treatment for 
j3 = 10, full quantum treatment for /3 = 10. In this example 
we used A'' = 18. 

Finally, by substituting Eqs. 1021) and (glj into Eq. 
we obtain 

'%^-5,jl-xn^~xC,tQ- (44) 

Short time results were obtained by integrating Eqs. 
%V,\\ . (|34|1 and 1)44(1 numerically using a fourth order 
Runge-Kutta procedure. Those numerical simulations 
should reproduce the degenerate model in the limit (3 <^ 
1. Note however that the short-time approximation 138|) 
is slightly different from that made in order to obtain 
Ea. (|19|l . In particular, Eq. (|38|l ignores a contribution 
linear in fit and as a consequence {fii,) increases more 
slowly in the approximate non-degenerate simulation so 
that for small (3 it agrees only within ~ 10% with the de- 
generate model once 5% of the initial population of atom 
pairs has been converted into molecules. This is illus- 
trated in Fig. which plots the average molecule number 
as a function of time. 

Although calculations with larger particle number are 
computationally tractable for short times, we have re- 
stricted our calculations to 18 atom pairs in order to al- 
low an exact numerical treatment of the full quantum 



non-degenerate model. Figure shows from left to right 
the degenerate analytical solution H21|) , the full quantum 
result with (3 = 0.1, the short time result with (3 = 0.1, 
short time result with /3 = 10 and the full quantum re- 
sult for f3 = 10, respectively. The full quantum result 
with (3 = 0.1 is indistinguishable from the full quantum 
treatment in the degenerate model, as expected. These 
simulations illustrate in particular that the formation of 
molecules is suppressed for increased (3. This is a conse- 
quence of the fact that a significant fraction of the atom 
pairs are detuned from resonance and therefore converted 
more slowly and incompletely into molecules. 

V. SUMMARY AND OUTLOOK 

We have shown that the early stages of molecular 
dimer formation from fermionic atoms are characterized 
by large fluctuations that reflect the quantum fluctua- 
tions in the initial atomic state. In contrast, the re- 
verse process of dissociation of a condensate of molec- 
ular dimers is largely deterministic. The reason for this 
asymmetry can be traced to the fact that in contrast to 
a quantum-degenerate fermionic system, the initial state 
of the molecular condensate is well described by a mean- 
field theory, that is, it is largely classical and relatively 
devoid of quantum fluctuations. 

As long as the atom-molecule coupling is dominant, 
f3 ^ 1, the kinetic energies are unimportant and a de- 
generate model can accurately describe the molecule for- 
mation. This is confirmed by our numerical simulations 
which show that the results of the degenerate model and 
the full quantum results are indistinguishable for /3 < 0.1. 
Small deviations appear for /3 w 1, and the creation of 
molecules is dramatically suppressed for higher ratios. 

Future work will extend these considerations to a more 
realistic multimodc description of the association process 
as well as a more detailed description of the two-body 
physics, which should in particular include the dressing 
of molecules by the open channel atomic pairs that is 
important in Feshbach sweeps. We will also extend our 
model to the case of fermionic molecules. 

This work is supported in part by the US Office of 
Naval Research, the NSF, the US Army Research Office, 
NASA, and the Joint Services Optics Program. 
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